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SYSTEM, METHOD, AND COMPUTER PROGRAM PRODUCT 

FOR SMOOTHING 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 [0001] This application claims the benefit of the filing 

date of United States Provisional Patent Application 

60/498,085 filed August 26, 2003, which is hereby 
incorporated by reference. 



10 TECHNICAL FIELD OF THE INVENTION 

[0002] The present invention is directed, in general, to 
finite element processing. 

BACKGROUND OF THE INVENTION 

15 [0003] With a very strong legacy of finite element data 

models, in most engineering industries, more and more 
legacy FEM data is being read in as building blocks for new 
designs. This is fast changing the engineering design world 
to an analysis-driven design as opposed to traditional CAD- 

20 based design. 

[0004] In this light, one of the biggest challenges that 
the designer faces today is to create geometry from legacy 
FEM data that is partially altered. The first step towards 
building surfaces or geometry abstractions from the input 
25 mesh is to smooth shell meshes in 3D space. Very few 3D 

smoothers, known today, can meet this challenge. Such a 
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smoother could also be looked at as a preliminary "mesh 
morphing" technique . 

[0005] The simplest mesh smoothing technique is 
Laplacian smoothing where a node is moved to the centroid 
5 of its neighboring vertices. This method operates 

heuristically, and has no control on mesh quality and often 
throws nodes outside concave domains producing inverted or 
invalid elements. 

[0006] Several variants of the Laplacian smoothing 
10 technique include "smart" Laplacian methods where element 

distortion metrices are checked before the node is moved. 
Length or area weighted Laplacian methods are able to 
influence repositioning with initial location. In concave 
domains, however, these techniques still fail to produce 
15 stable meshes. Sometimes, area weighted methods help sense 

inverted elements and can be leveraged by smart Laplacian 
approaches to prevent moving such nodes. Haber et al . 
proposed a family of Isoparametric-Laplacian mesh control 
techniques, meant specially for transfinite meshes. These 
2 0 schemes tend to produce better skewed quadrilaterals but 

have most of the other Laplacian disadvantages. 

[0007] In optimization based smoothing, nodes are not 
moved based on a heuristic algorithm such as Laplacian 
smoothing, they are moved to minimize a given distortion 
25 metric. The distortion metrics are related to the max/min 

included angles of elements, element skew value, element 
aspect ratio, element area, element edge length, etc. 

[0008] One of the early optimization techniques was 
developed by Parthasarathy and Kodiyalam. They solved a 



- 2 - 



DOCKET NO. 05-03-002 



PATENT 



non-linear optimization problem in an effort to repair 
quad-tree and octree meshes. Shephard and Georges reported 
similar findings. Amenta et al used linear programming 
techniques to solve triangular meshes locally. Jacquotte 
5 and Coussement developed an optimization based approach for 

both 2D and 3D structural grids. Freitag et al proposed a 
local optimization technique for 2D triangular meshes that 
can serve as the core of an efficient parallel algorithm. 
Later, Freitag and Oliver-Gooch extended that to 3D grids. 

10 These methods produced extremely good quality meshes but 

involved element repair work like edge -swapping, bad 
element collapse etc. Moreover the efficiency of these 
optimization based smoothers are about five to 40 times 
slower than Laplacian smoothing. Both Freitag and Canann et 

15 al later combined several smart Laplacian methods with 

optimization-based techniques to create hybrid algorithms 
to improve efficiency. These methods, however, fail to 
recognize and preserve mapped meshes. 

[0009] Non-Laplacian, physics-based or non- iterative, 
2 0 direct smoothing algorithms have also emerged in the recent 

past. Lohner et al used a spring mounted system between 
nodes to smooth. Shimada, in his bubble packing model, 
proposed a method where close nodes repel and distant ones 
attract each other. Bossen and Heckbert proposed an 
25 exponential function with similar properties. These models 

work in 2D space and are not general purpose. A recent 
direct, non- iterative approach computes artificial 
stiffness matrices for the mesh to smooth and tries to 
minimize the strain-energy of the system. While the results 
30 are interesting, such methods tend to loose their 

computational advantage for large models. Tarn and Armstrong 
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developed an integer programming based mesh control 
algorithm for mapped quadrilateral and hexahedral meshes.. 
This method works on a collection of connected subregions 
and is quite special purpose. Zhou and Shimada have 
5 recently proposed an angle smoother in 2D that tends to 

mount torsion springs between nodes and minimize the system 
energy. Results prove its merit on certain concave domains 
and mixed meshes. Winslow smoothing is another efficient 
technique to reposition nodes of predominantly structured 
10 elliptic meshes. Knupp ■ s investigation in this area 

provides important results. Knupp also proposes "condition 
numbers" as a good yardstick for measuring mesh quality. 

[0010] There is a need in the art for a system, process, 
and computer program product for improved smoothing. 
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SUMMARY OF THE INVENTION 

[0011] The preferred embodiment provides a system, 
method, and computer program product for a hybrid, 
variational, user-controlled, 3D mesh smoothing for 
5 orphaned shell meshes. The smoothing model is based on a 

variational combination of energy and equi-potential 
minimization theories, A variety of smoothing techniques 
for predicting a new location for the node -to- smooth are 
employed. Each node is moved according to a specific 

10 smoothing algorithm so as to keep element included angles, 

skew and distortion to a minimum. The variational smoother 
selection logic is based on nodal valency and element 
connectivity pattern of the node to smooth. Results show 
its consistency with both quadrilateral and quad-dominant 

15 meshes with a significant gain over conventional Laplacian 

schemes in terms of mesh quality, stability, user control 
and flexibility. 

[0012] The foregoing has outlined rather broadly the 
features and technical advantages of the present invention 

20 so that those skilled in the art may better understand the 

detailed description of the invention that follows. 
Additional features and advantages of the invention will be 
described hereinafter that form the subject of the claims 
of the invention. Those skilled in the art will appreciate 

25 that they may readily use the conception and the specific 

embodiment disclosed as a basis for modifying or designing 
other structures for carrying out the same purposes of the 
present invention. Those skilled in the art will also 
realize that such equivalent constructions do not depart 
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from the spirit and scope of the invention in its broadest 
form. 

[0013] Before undertaking the DETAILED DESCRIPTION OF 
THE INVENTION below, it may be advantageous to set forth 
5 definitions of certain words or phrases used throughout 

this patent document: the terms "include" and "comprise," 
as well as derivatives thereof, mean inclusion without 
limitation; the term "or" is inclusive, meaning and/or; the 
phrases "associated with" and "associated therewith," as 

10 well as derivatives thereof, may mean to include, be 

included within, interconnect with, contain, be contained 
within, connect to or with, couple to or with, be 
communicable with, cooperate with, interleave, juxtapose, 
be proximate to, be bound to or with, have, have a property 

15 of, or the like; and the term "controller" means any 

device, system or part thereof that controls at least one 
operation, whether such a device is implemented in 
hardware, firmware, software or some combination of at 
least two of the same. It should be noted that the 

20 functionality associated with any particular controller may 

be centralized or distributed, whether locally or remotely. 
Definitions for certain words and phrases are provided 
throughout this patent document, and those of ordinary 
skill in the art will understand that such definitions 

25 apply in many, if not most, instances to prior as well as 

future uses of such defined words and phrases. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] For a more complete "understanding of the present 
invention, and the advantages thereof, reference is now 
made to the following descriptions taken in conjunction 
5 with the accompanying drawings, wherein like numbers 

designate like objects, and in which: 

[0015] Figure 1 depicts a typical n-quad grid 
illustrating a n-quadrilateral connection in accordance 
with an embodiment of the present invention; 

10 [0016] Figures 2A-2C illustrate some of the aspects of 

an isoparametric-laplacian smoothing process in accordance 
with an embodiment of the present invention; 

[0017] Figure 3 shows a spring-mounted system in 
accordance with an embodiment of the present invention; 

15 [0018] Figures 4A-4C illustrate the effectiveness of the 

disclosed system and method in improving element included 
angles; 

[0019] Figure 5 depicts a node connected to 4 
quadrilateral elements ; 

20 [0020] Figures 6A and 6B depict a merged single quad in 

accordance with a preferred embodiment; 

[0021] Figures 7A and 7B depict a single quad after 
constrained smoothing in accordance with a preferred 
embodiment ; 
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[0022] Figure 8 shows a 2D quad-grid where the node to 
smooth is constrained to move within the bounding box of 
the three centroids of the three quads, in accordance with 
an embodiment of the present invention; 

5 [0023] Figures 9A-9C illustrate some benefits of the 

disclosed variational smoother; 

[0024] Figures 10A-10C show a 5X6 mapped all quad mesh 
in a concave region, in accordance with an embodiment of 
the present invention; 

10 [0025] Figure 11 shows a portion of an exemplary car 

body panel that has an initial distorted mesh; 

[0026] Figure 12 shows a mesh smoothed by a variational 
smoothing process in accordance with a preferred 
embodiment ; 

15 [0027] Figure 13 depicts a block diagram of a data 

processing system in which a preferred embodiment can be 
implemented; and 

[0028] Figure 14 depicts a flowchart of a process in 
accordance with a preferred embodiment. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0029] FIGURES 1 through 12, discussed herein, and the 
various embodiments used to describe the principles of the 
present invention in this patent document are by way of 
5 illustration only and should not be construed in any way to 

limit the scope of the invention. Those skilled in the art 
will understand that the principles of the present 
invention may be implemented in any suitably arranged 
device. The numerous innovative teachings of the present 
10 application will be described with particular reference to 

the presently preferred embodiment. 

[0030] The preferred embodiment provides a system and 
method for a hybrid, variational, user-controlled, 3D mesh 
smoothing algorithm for orphaned shell meshes. The 

15 smoothing model is based on a variational combination of 

energy and equi-potential minimization theories. A variety 
of smoothing techniques for predicting a new location for 
the node -to- smooth are employed. Each node is moved 
according to a specific smoothing algorithm so as to keep 

20 element included angles, skew and distortion to a minimum. 

The variational smoother selection logic is based on nodal 
valency and element connectivity pattern of the node to 
smooth. Results show its consistency with both 
quadrilateral and quad-dominant meshes with a significant 

25 gain over conventional Laplacian schemes in terms of mesh 

quality, stability, user control and flexibility. 

[0031] The preferred embodiment discloses a new 
variational shell mesh smoothing technique in 3D space that 
symbiotically combines old and some new methods. The 
3 0 variational algorithm smoothes each node according to a 
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specific smoothing technique. The smoothing method 
selection depends on nodal valency and connectivity pattern 
of the concerned node and is largely determined based on 
different smoothing techniques on a variety of mesh 
5 patterns. The algorithm used in various embodiments is 

iterative; works on ID, 2D and 3D space; is almost as 
efficient as Laplacian smoothing; gives several controls to 
the user; does not destroy mapped/structured meshes or mesh 
regions; does not produce inverted elements; does not move 
10 boundary nodes, or nodes on hard points or nodes with 

loads; and improves element included angles, average 
element skew and hence mesh quality. 

[0032] Many users handling geometry-associated meshes 
desire to move nodes off the surface in an effort to reduce 
15 element warp and skew. For geometry associated meshes, this 

algorithm allows the user to move nodes off the surface 
within a desired tolerance, if needed. 

[0033] The preferred embodiment uses different smoothing 
method for each node in the mesh, including an incenter- 

2 0 based approach for triangles, isoparametric -Laplace for 

general quad-only meshes, equipotential (Winslow) smoothing 
for. mapped regions, and combined incenter and Laplacian 
smoothing for free mixed meshes. Further, the process does 
a quick region check for each node, optionally performs a 

25 smart interior angle screening, and optionally constrains 

the node to move within a specified tolerance. 
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[0034] The governing equation according to one 

preferred embodiment, for repositioniong node i connected 
to N elements, is 

N 

Pi=ZFn(C,V)*On(C,V) (1) 
n= 1 

where Pi' = New position of node i, Fn = Variational weight 
factor for n-th element an = Positional function for n-th 
element, C denotes the connectivity pattern of the node, 
while V indicates its valency. 

[0035] The hybrid smoother uses various repositioning 
schemes for different types of mesh units. These schemes 
are disussed below. 

[0036] Conventional Laplacian methods tend to ignore the 
effect of neighboring nodes that are not directly connected 
to the node -to -smooth. Furthermore, being heuristic and 
completely independent of element topology, these 
techniques do not guarantee low skewness for quad elements. 
This is a strong requirement for most structural 
engineering applications. Isoparametric mesh generation 
fundamentals, derived from the shape functions of 
quadrilateral isoparametric elements, predict location of 
nodes connected to quadrilateral elements. The preferred 
embodiment makes a concious effort to generate grids that 
adequately represent bounding contours and tend to produce 
squarish quad elements that have low skew. When this 
technique is combined with Laplacian grid generation 
procedures, a more uniform mesh results that have a lower 
skew value than ordinary Laplacian grids. Such an 
Isoparametric -Laplacian method is used here as a mesh 
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smoothing technique for n-quad connections where n ! = 2 or 
4. A typical n-quad grid is shown in Figure 1, which 
illustrates a n-quadrilateral connection (n != 2,4). i-j-k- 
1 denotes an average quad, i is the node to be smoothed. 

[0037] From the origin of the isoparametric co-ordinate 
systems for quadrilateral elements, the new location of 
node i as shown in Figure 1 can be expressed as 

l N 

Pi' = x wn (Pnj + Prit - rPn k ) (2) 

N(2-w) n=l 

where the variational weighting factor and the positional 
function for each element can be expressed as 

Wn 

Fn = ; Qn = (Pnj + Pnj - rPn k ) (2a) 

N(2-r) 

and N is the number of elements connected to node I, i-j-k- 
1 represents an average connected quad, r is the coupling 
factor between Laplacian and isoparametric methods, Pnj, 
Pnk, Pnl, represent the position vectors of the j-th, k-th 
and 1-th nodes of the n-the connected quad respectively, 
Pi' represents the new location of the node to smooth, and 
Wn represents weight factors for each connecting element n 
such that 

N 

I Wn = 1.0; 
n=0 

[0038] The weight factors can be constructed in numerous 
possible ways. Some popular examples are length-weighing, 
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area-weighing, weighing elements based on the included 
angle they make at the node concerned, etc. 

[0039] Note that when r = 0, equation (1) reduces to 
Laplace smoothing. When r = 1.0, a pure isoparametric grid 
is produced with quad elements showing very low skewness, 
but the nodal lines of the mesh become zig-zag. Experience 
with scheme has proven that r = 0.5 results in good quality 
meshes with an overall skewness that is almost invariably 
better than the Laplacian variants. 

[0040] Figures 2A-2C illustrate some of the aspects of 
the isoparametric- laplacian smoothing scheme discussed 
above. As is evident from Figure 2A, the simple Laplacian 
smoother does just a fairly good job smoothing the 5-valent 
quad connection but the iso-laplace method produces better 
quad included angles. The isoparametric smoother, on the 
contrary, tends to produce a lot more skew elements while 
trying preserve quad-shapes. 

[0041] Figure 2A illustrates Laplace smoothing, Figure 
2B illustrates isoparametric smoothing, and Figure 2C 
illustrates Iso-Laplace smoothing. Each of these shows a 
dark smoothed mesh overlaying a thinner original mesh. 

[0042] A new Incenter-Based Smoother scheme is used for 
triangles only. Figure 3 shows a spring-mounted system 
where the node -to- smooth is assumed to be connected by 
springs to the incenters of the connected triangles. An 
energy minimization similar to Laplace smoothing leads to 
the following expressions. 
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[0043] The new position of node i is obtained as 

N 

Pi' = Pi + I W n (P n - PJ (3) 

n=l 

5 

where Pi (x, y, z) is the position vector of node I, Pn(x, 
y, z) is the incenter vector of element n, and N represents 
the total number of elements connected to node i . The 
weighing factors Wn account for the initial stiffness of 
10 the springs (i.e the effect of initial position of node i) . 

Equation (2) is ambivalent with a length weighted Laplacian 
method. The weight factors, in this case, can be computed 
from the length of the incenter vectors as 

llv.ll 

15 W n = (3a) 

N 

Z||v n || 
n=l 

20 where 

v n = P.- Pi (3b) 

Vn is the distance from the concerned node to the incenter 
of the n-th element. 

[0044] Figure 3 shows an all -triangle connection. 
25 Springs are connected from the node -to- smooth (i) to the 

Incenters of the connected triangles. 

[0045] Figures 4A-4C illustrate the effectiveness of the 
disclosed system and method in improving element included 
angles. Smoothing results are shown for a unit mesh 
30 comprising three triangles connected to a single node. 
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[0046] Figure 4A shows three triangles connected to one 
interior node. Figure 4B shows the same after Laplace 
smoothing. Figure 4C shows the same after incenter 
smoothing. The mininum and maximum angles are compared in 
5 Table 1. 



Table 1. Comparison of element included angles 



Smoothing Method 


Minimun Element Angle 


Maximum Element Angle 


Original Mesh 


3.7 


171.2 


Laplacian smoothing 


17.3 


137.6 


Angle Smoothing [201 


16.9 


133.2 


Disclosed variational smoothing 


17.4 


125.7 



10 [0047] The equipotential smoothing scheme applies to 

nodes that are connected to n quads (where n = 4) or has 8 
neighboring nodes. Figure 5 depicts a node connected to 4 
quadrilateral elements. Neighboring nodes 1,2.. 8 are shown. 

[0048] The governing equation for equipotential 
15 (Winslow) smoothing can written for node i as 

aPi u - 2pPi^ n + yPi nn = 0; (4) 

where x,h are logical variables that are harmonic in 
nature, while a, b, g are constant coefficients that 
depend on the problem. 

20 [0049] Figure 5 shows a four-quad grid connected to node 

i. From equations 1 & 4, after some simplification, one 
obtains an equation similar to equation 3 that describes 
the new location of the node to smooth, where N = 8. The 
weighing factors of the 8 neighboring nodes are given by 
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W,=-p/2,W 2 =a,W 3 =p/2, W 4 =y,W 5 =-p/2 J W 6 =a , W 7 =p/2,W 8 =y (4a) 

where 

a = x p 2 + y p 2 + Zp 2 (4b)P = Xp x q + y p y q + (4c) 
Y = x q 2 + y q 2 + z q 2 (4d)and 

x p = (x 2 -x 6 )/2, y p =(y 2 -y 6 )/2, Zp = (z 2 -Zfi)/2 (4e) 
x q = (x 8 -x 4 )/2, y q =(y 8 -y 4 )/2, Zq = (z8-z 4 )/2 (4f) 

[0050] A Combined Incenter-IsoLaplacian smoothing for 

mixed meshes scheme applies to nodes that are connected to 
both quads and triangles and do not satisfy the Winslow 
criteria. For such connections, a Isoparametric-Laplacian 
approach is used for quads and the incenter-based smoothing 
approach is used for triangles. The region check is turned 
on to prevent mesh distortion in concave domains. 

[0051] A quad meshes with bivalent nodes scheme is used 
when we come across a node with a rare valency of 2 

(connected to 2 quads) . One of these quads is invariably 
an arrowhead shaped element which is impossible to cure by 
any smoothing method. For such bivalent nodes, the node to 
smooth is deleted and the two connecting quads are merged 
to form a single quad as shown in Figures 6A & 6B. Figure 
6A shows elements with bi-valent node, and Figure 6B shows 
the Bivalent-node removed, where elements merge. 

[0052] Angle check: During smoothing, element angles 
connected to the concerned node are checked to see if they 
are within the angle limits, or the overall skewness of the 
connected elements are checked. If all elements connected 
to the node pass the user-defined criteria, the concerned 
node is not smoothed. This feature is optional, and when 
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turned on, the smoother works like a smart Laplacian. Node 
i, is thus smoothed, only if any included angle of any 
connected element fails the user-set angular limits at this 
node. This condition can be mathematically expressed as 

a (Pi) ; if Omax < | | j . i, N < Gmin 00 (5) 

[0053] where a is the smoothing operator, a (Pi) 
indicates, node i is smoothed, N denotes number of elements 
connected to node I, 0max, Gmin are the element allowable 
angular limits, and otji denotes the included angle of 
element element j at node i. 

[0054] Constrained Movement: Since this smoother is most 
beneficial to users who are trying to modify existing 
meshes, add features and build geometry from them, many 
users would prefer the to move nodes in a constrained 
manner. For example, they may not want their mesh nodes to 
move off their current location by more than a certain 
amount. Once the user specifies such a node movement 
tolerance dTol, the nodes are moved from their present 
location such that they always reside within the sphere of 
radius R = dTol as shown in Figures 7A and 7B. 

[0055] The new location of node i (Pi") is given by 
Pr = rP k + (!-r)Pi ; for r < 1 (6a) 
P r = rP k ;forr>l (6b) 

where 



dTol 
r = 



(6c) 
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and Pi, Pi", and Pk are the node locations before 
smoothing, after constrained smoothing and after 
unconstrained smoothing respectively. 

[0056] Region Check: While the hybrid spring system 
tries to reposition the node so as to relax the worst 
angles of the connected elements, in certain situations, it 
is still bound to suffer from the so-called "Laplacian 
distortion" effect, where, interior nodes might move 
outside concave boundaries. To ensure that the moved node 
lies inside the domain of its neighbors, a region, is 
defined by a bounding box of the neighbors. The node to 
smooth should not be moved outside this bounding box. 

[0057] Let the bounding box of the neighbors be 
expressed as 

fi(X 1 ,X 2 ,Y 1 ,Y 2) Z 1> Z 2 ) (7a) 

[0058] where XI = Xmin of all the initial X's of the 
centroid (for quads) /incenters (for triangles) of the 
neighboring elements, X2 = Xmax of all the initial X's of 
the centroid (for quads) /incenters (for triangles) of the 
neighbors, and Yl, Y2 , Zl and Z2 are expressed identically. 

[0059] For the node movement to be valid the following 
conditions need to be met 

[0060] New node location should lie inside bounding box 
n(Xl, X2, Yl, Y2, Zl, Z2) 

X2 > Px > XI 
Y2 > Py > Yl 

Z2 > Pz > Zl (7b) 
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where Px,Py,Pz denote the new coordinates of the node. 

[0061] Figure 8 shows a 2D quad-grid where the node to 
smooth is constrained to move within the bounding box of 
the three centroids of the three quads. 

5 [0062] Convergence: Since most of the methods used in 

the variational smoothing algorithm are heuristic in 
nature, the error at each step can be defined by the root 
mean square of the positional disturbance of the mesh 
nodes. This can be expressed from equation (1) as 

10 Bj = Z n (P i .-P i ) 2 /n 

i = 0 

for the j-th iteration of a mesh that has n nodes. The 
smoother is assumed to have converged when the following 
criteria is met: 8j - Sj_i < O.OlSi. 

15 [0063] The variational logic is controlled by two 

factors: the valency of the node to smooth, and its element 
connectivity pattern. Table 2 can best express the rule set 
adopted. 

Table 2: Variational smoothing logic 

2 0 



Valency of node to 
smooth 


No of Quads 


No of Trias 


Smoothing Method 


2 


2 




Scheme E: node 
removed, 2 quads 
are combined into 
one 


3 


3 




SchemeA: 
Isoparametric- 
Laplace 


4 


4 




Scheme C: 
Equipotential 


n; n>4 


N 




Scheme A: 
Isoparametric- 
Laplace 
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n; n > 2 
n!=8 




n 


Scheme B* 
Incenter-based 


8 




8 


Scheme C: 
Equipotential 


n;n>2 
n = k + l 


K 


1 


Scheme D: 
Incenter-IsoLaplace 



[0064] Several unitary and large size shell meshes are 
smoothed using the disclosed model. Results are compared 
with simple Laplace or smart -Laplace smoothing cases or 
5 other existing smoothing algorithms. The ability of the 

smoother to clean grossly distorted meshes, maintain 
stability in concave and triangular regions, preserve 
mapped or structured meshes, improve interior angles and 
skewness, and converge efficiently is studied in the 
10 following subsections . 

[0065] Ability to preserve mapped -me she s : Most smoothing 
algorithms do not recognize mapped/structured grids. The 
angle smoother employs a Winslow smoothing algorithm [12] 
for nodes that are connected to 4 quadrilaterals or 8 
15 triangles. The benefits of the variational smoother are 

illustrated in Figures 9A-9C. Figure 9A shows a 4X6 mapped 
mesh on a concave region. 

[0066] Figure 9B shows the same mesh after Laplace 
smoothing. The structured nature of the mesh is lost. Edge 
20 elements have almost collapsed. 

[0067] Figure 9C shows the result of variational 
smoothing. All nodes are equipotentially smoothed 

resulting in a structured or mapped mesh pattern. No 
element is collapsed. 
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[0068] Ability to Work On Concave Boundaries: On 

concave domains Laplacian smoothing often pushes points 
outside the mesh boundary. The disclosed angle smoother 
prevents that either by employing a Winslow in such regions 
5 if the mesh qualifies for it or by a region check. Figures 

9A-9C show a 5X6 mapped mesh distorted by Laplacian 
smoothing but completely preserved by the current method. 

[0069] Distortions caused by Laplacian smoothing in 
concave domains by throwing nodes outside the mesh region 
10 and creating inverting elements can also be prevented by 

the region check feature of the variational smoother. 

[0070] Figure 10A shows a 5X6 mapped all quad mesh in a 

concave region. Figure 10B shows the result of Laplace 
smoothing, where nodes move outside the boundary resulting 
15 in inverted elements. Figure 10C shows the result of the 

di sclosed variational smoothing, where no nodes go outside 
the boundary and the mesh has not been significantly moved. 

[0071] Ability to Smooth Quad-dominant meshes: Figure 

11 shows a portion of a car body panel that has an initial 

20 distorted mesh which is cleaned by the variational smoother 

as shown in Figure 12 . Table 3 compares the performance of 
the angle smoother to smart -Laplacian smoothing. Once 
again, it is observed that the number of elements failing 
included angle check for the angle smoother is half of that 

25 produced by smart -Laplacian smoothing. 

[0072] Figure 11 shows a portion of a carbody panel with 
quad dominant mesh, partially mapped mesh with high-skew 
elements . 
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[0073] Figure 12 shows a mesh smoothed by the new 
variational smoother. All elements are smoothed. Resultant 
mesh shows mapped grids have been preserved, most high skew 
elements cured. 

5 Table 3. Comparison of mesh quality for a car-body panel 



Smoothing Method 


No. of elm 


No. of elem failing 


No. of elem failing 


No. of elem failing 




failing warp 


minimum Jacobian 


skew check 


included 


angle 




check 


check 




check 


Original mesh 


37 


89 


241 


156 


After Laplace smoothing 


5 


2 


89 


31 


After variational smoothing 


1 


2 


63 


15 



[0074] Laplacian smoothing or its variants have certain 
inherent shortcomings that cannot be cured by a single 

10 alternate smoothing algorithm. More over, most schemes are 

restricted to either two-dimensional or monolithic grids. 
The variational smoother of the preferred embodiment uses a 
hybrid smoothing mechanism where different smoothing (both 
conventional and innovative) algorithms are used for each 

15 node depending on its connectivity. For all triangular 

meshes, an incenter-based method is disclosed that performs 
better and less destructively than the Laplacian method. 
For all quadrilateral meshes, an Isoparametric -Laplace 
method is used that leads to faster convergence and lower 

20 skewness. Winslow smoothing is employed for mapped-grid 

units that preserves structured grids and is non- 
destructive in concave domains. The method is made more 
flexible by adding an user control on the element check 
criteria (included angle or skew) and is further reinforced 

25 by a quick region check that often ensures valid meshes in 

skew and concave contours. Several results exemplify the 
performance of the disclosed variational smoother in 
comparison to a conventional Laplacian, smart -Laplacian and 
other schemes . 



- 22 - 



DOCKET NO. 05-03-002 



PATENT 



[0075] Figure 13 depicts a block diagram of a data 
processing system in which a preferred embodiment can be 
implemented. The data processing system depicted includes a 
processor 1302 connected to a level two cache/bridge 1304, 
which is connected in turn to a local system bus 1306. 
Local system bus 1306 may be, for example, a peripheral 
component interconnect (PCI) architecture bus. Also 
connected to local system bus in the depicted example are a 
main memory 1308 and a graphics adapter 1310. 

[0076] Other peripherals, such as local area network 
(LAN) / Wide Area Network / Wireless (e.g. WiFi) adapter 
1312, may also be connected to local system bus 1306. 
Expansion bus interface 1314 connects local system bus 1306 
to input/output (I/O) bus 1316. I/O bus 1316 is connected 
to keyboard/mouse adapter 1318, disk controller 1320, and 
I/O adapter 1322. 

[0077] Also connected to I/O bus 1316 in the example 
shown is audio adapter 1324, to which speakers (not shown) 
may be connected for playing sounds. Keyboard/mouse adapter 
1318 provides a connection for a pointing device (not 
shown), such as a mouse, trackball, trackpointer , etc. 

[0078] Those of ordinary skill in the art will 
appreciate that the hardware depicted in Figure 13 may vary 
for particular. For example, other peripheral devices, such 
as an optical disk drive and the like, also may be used in 
addition or in place of the hardware depicted. The depicted 
example is provided for the purpose of explanation only and 
is not meant to imply architectural limitations with 
respect to the present invention. 



- 23 - 



DOCKET NO. 05-03-002 



PATENT 



[0079] A data processing system in accordance with a 
preferred embodiment of the present invention includes an 
operating system, optionally employing a graphical user 
interface. When a graphical user interface is employed, the 
5 operating system permits multiple display windows to be 

presented in the graphical user interface simultaneously, 
with each display window providing an interface to a 
different application or to a different instance of the 
same application. A cursor in the graphical user interface 
10 may be manipulated by a user through the pointing device. 

The position of the cursor may be changed and/or an event, 
such as clicking a mouse button, generated to actuate a 
desired response. 

[0080] One of various commercial operating systems, such 
15 as Unix, Linux, or a version of Microsoft Windows™, a 

product of Microsoft Corporation located in Redmond, Wash, 
may be employed if suitably modified. The operating system 
is modified or created in accordance with the present 
invention as described. 

20 [0081] Figure 14 depicts a flowchart of a process in 

accordance with a preferred embodiment. Here, a finite 
element model is first loaded into a data processing system 
(step 1405) . It should be noted that the model is not 
necessarily loaded as a graphic image, but can be loaded as 

25 data that describes the model, as known to those of skill 

in the art. Next, a node of the model is selected for 
smoothing (step 1410) . Further, it should be noted that 
while the preferred embodiments are drawn to finite element 
models, these methods can also be applied to more simple 

30 graphic models. 
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[0082] The selected node is then analyzed. The nodal 
valency of the node is determined (step 1415) , and the 
element connectivity pattern of the node is determined 
(step 1420) . 

[0083] Next, the node is smoothed according to the nodal 
valency and element connectivity pattern. If the element 
connectivity pattern is a triangle, incenter-based 
smoothing is performed (step 1425) ; if the element 
connectivity pattern is quad-only mesh, isoparametric- 
Laplace smoothing is performed (step 143 0) ; if the element 
connectivity pattern is a mapped region, equipotential 
smoothing is performed (step 1435) ; and if the element 
connectivity pattern is a free-mixed mesh, combined 
incenter and Laplacian smoothing is used (step 1440) . 

[0084] Preferably, the smoothing of the node is 
performed according using: 

N 

Pi' = lFn(C,V)*On (C,V) 
n= 1 

where i is the node to be smoothed, i is connected to N 
elements, Pi' is the new position of node i, Fn is the 
variational weight factor for n-th element Qn is the 
positional function for n-th element, C denotes the 
connectivity pattern of the node, and V indicates the 
valency of the node . 

[0085] Next, the system optionally performs an interior 
angle screening (step 1445) , according to user-specified 
criteria, and optionally contrains the node to move within 
a predetermined tolerance (step 1450) . 
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[0086] The modified model is then stored (step 1455) . 

[0087] Further details of mesh smoothing and processing 
techniques can be found in the following, all of which are 
hereby incorporated by reference: Amenta, N., .Bern, M. and 
5 Eppstein, D. (1997) Optimal point placement for mesh 

smoothing. Proc . 8th ACM-SIAM Symp. On Discrete Algorithms, 
New Orleans, LA; Balendra, B. (1999) A direct smoothing 
method for surface meshes. Proc. 8th Intl. Meshing 
Roundtable, Sandia National Laboratories, 303-308; Bossen, 

10 F.J. and Heckbert, P.S. (1996) A pliant method for 

anisotropic mesh generation. Proc. 5th Intl. Meshing 
Roundtable, Sandia National Laboratories, 63-74; Canann, 
S.A., Tristano, J.R. and Staten,M.L. (1998) An approach to 
combined Laplacian and Optimization-based smoothing for 

15 triangular, quadrilateral and quad-dominant meshes. 7th 

Intl. Meshing Rountdable, Sandia National 

Laboratories; Field, D. (1988) Laplacian smoothing and 
Delauney triangulations . Comm. App .Num.Meth, 4, 709-712; 
Freitag, L., Jones, M. and Plassmann, P. (1995) An 

20 efficient parallel algorithm for mesh smoothing. Proc. 4th 

Intl. Meshing Roundtable, Sandia National Laboratories, 47- 
58; Freitag, L. and Oliver-Gooch, C. (1997) A comparison of 
tetrahedral mesh improvement techniques. Proc. 6th Intl. 
Meshing Roundtable, Sandia National Laboratories, 24 9-2 62; 

25 Freitag, L. (1997) On combining Laplacian and optimization- 

based mesh smoothing techniques. ASME, 220, 37-43; Haber, 
R. Shephard, M. S., Abel, J. F. , Gallagher, R.H., and 
.Greenberg, D.P. (1981) A general two-dimensional, 
graphical finite element preprocessor utilizing discrete 

30 transfinite mappings. IJNME 17, 1015-1044; Ives,D. (2000) 

Unstructured boundary layer grid generation. Proc. 7th 
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Intl.Conf. Num. Grid Generation in Computational Field 
Simulations, British Columbia, Canada, 13-19; Jacquotte, 
O.P. and Coussement, G. (1992) Structured mesh adaptation: 
space accuracy and interpolation methods. Comp. Meth. App. 
5 Mech. Engg, 101 397-432; Jones, R.E. (1977) A self- 

organizing mesh generation program. J. Pressure Vessel Tech, 
ASME, 96 (3), 193-199; Lee, C.K. and Lo, S.H. (1994) A new 
scheme for the generation of a graded quadrilateral mesh. 
Comp. Struct. 52(5), 847-857; Lo,S. (1985) A new mesh 

10 generation scheme for arbitrary planar domains. IJNME, 21, 

1403-1426; Lohner, R. , Morgan, K. and Zienkiewicz , 0 . C . 
(1986) Adaptive grid refinement for the compressible Euler 
equations. Accuracy estimates and adaptive refinements in 
finite element computations, Ed. I. Babuska, 

15 O.C. Zienkiewicz, J.Gago and E.R.de A. Oliviera, Wiley, 281- 

297; Parthasarathy, V. and Kodiyalam, S. (1991) A 
constrained optimization approach to finite element mesh 
smoothing. J. Fin. Elm. Anal .Des, 9, 309-320; Shephard, M.S. 
and Georges, M. (1991) Automatic three-dimensional mesh 

20 generation by the finite octree technique. IJNME, 32, 709- 

749: Shimada,K. (1993) Physically-based mesh 

generation: automated triangulation of surfaces and volumes 
via bubble packing. PhD Thesis, ME Dept., MIT; Tarn, T.K.H. 
and Armstrong, C.G. (1993) Finite element mesh control by 

25 integer programming. IJNME, 36, 2581-2606; Zhou, T. and 

Shimada,K. (2000) An angle-based approach to two-dimensional 
mesh smoothing. Proc . 9th Intl. Meshing Roundtable, Sandia 
National Laboratories, 373-384; and Knupp, P. (1998) 
Winslow Smoothing On Two-Dimensional Unstructured Meshes. 

30 Proc. 7th Intl. Meshing Roundtable, Sandia National 

Laboratories, 449-457; Knupp, P. (1999) Applications of 
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Mesh Smoothing: Copy, Morph, and Sweep on Unstructured 
Quadrilateral Meshes", IJNME, Wiley, Vol 45, 37-45. 

[0088] Those skilled in the art will recognize that, for 
simplicity and clarity, the full structure and operation of 
all data processing systems suitable for use with the 
present invention is not being depicted or described 
herein. Instead, only so much of a system and method as is 
unique to the present invention or necessary for an 
understanding of the present invention is depicted and 
described. The remainder of the construction and operation 
of the system and method may conform to any of the various 
current implementations and practices known in the art. 

[0089] It is important to note that while the present 
invention has been described in the context of a fully 
functional system, those skilled in the art will appreciate 
that at least portions of the mechanism of the present 
invention are capable of being distributed in the form of a 
instructions contained within a machine usable medium in 
any of a variety of forms, and that the present invention 
applies equally regardless of the particular type of 
instruction or signal bearing medium utilized to actually 
carry out the distribution. Examples of machine usable 
mediums include: nonvolatile, hard-coded type mediums such 
as read only memories (ROMs) or erasable, electrically 
programmable read only memories (EEPROMs) , user-recordable 
type mediums such as floppy disks, hard disk drives and 
compact disk read only memories (CD-ROMs) or digital 
versatile disks (DVDs) , and transmission type mediums such 
as digital and analog communication links. 
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[0090] Although an exemplary embodiment of the present 
invention has been described in detail, those skilled in 
the art will understand that various changes, 
substitutions, variations, and improvements of the 
invention disclosed herein may be made without departing 
from the spirit and scope of the invention in its broadest 
form. 

[0091] None of the description in the present 
application should be read as implying that any particular 
element, step, or function is an essential element which 
must be included in the claim scope: THE SCOPE OF PATENTED 
SUBJECT MATTER IS DEFINED ONLY BY THE ALLOWED CLAIMS. 
Moreover, none of these claims are intended to invoke 
paragraph six of 35 USC §112 unless the exact words "means 
for" are followed by a participle. 



